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Abstract 

A  mathematically  robust  method  is  derived  that  allows  the  SAR  of  an 
electro-optical  system  to  be  estimated  when  a  point  source  image  falls  onto 
a  sensor  array  of  discrete  picture  elements  (pixel).  The  derivation  is  based 
on  a  geometrical  analysis  of  image  and  sensor  element  configuration.  This 
method  allows  the  optimal  pixel  size  to  be  selected  to  maximize  the  expected 
SNR  for  any  point  source  Image  once  It's  diameter  is  known.  In  this  report 
the  method  is  used  to  determine  the  probability  of  an  image  falling  into  one, 
two,  .  .  .,  six  pixels  as  a  function  of  Image  diameter  and  pixel  size,  and 
the  corresponding  expected  maximum  signal.  Examples  are  given  using  this 
method  to  estimate  the  SNR  for  an  electro-optical  system. 
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I.  INTRODUCTION 

The  analysis  of  the  performance  of  an  electro-optical  system  to  detect 
point  source  Images  Is  greatly  simplified  by  the  assumption  that  the  Images 
each  fall  In  one  sensor  array  picture  element1  (pixel).  This  assumption  Is 
not  valid  when  hundreds  of  star  Images  and  one  or  more  unresolved  satellite 
Images  (due  to  reflected  sunlight)  fall  on  a  focal  plane  of  thousands  of 
pixels,  because  although  some  of  the  Images  will  fall  entirely  In  one  pixel 
many  will  fall  in  two  or  more  pixels.  When  an  Image  falls  In  one  pixel,  the 
maximum  electrical  signal  generated  Is  proportional  to  the  total  Image  photons, 
but  when  an  Image  falls  In  more  than  one  pixel  the  maximum  electrical  signal 
becomes  a  function  of  the  distribution  of  the  Image  photons  In  the  various 
pixels.  One  of  the  pixels  will  contain  most  of  the  Image  photons  and  will 
give  the  maximum  electrical  signal.  The  maximum  electrical  signal  is  of  para¬ 
mount  Importance  because  If  the  maximum  signal  crosses  a  threshold  the  Image 
can  be  detected,  but  If  the  maximum  signal  fails  to  cross  a  threshold  the 
image  will  not  be  detected.  The  maximum  signal  resulting  from  an  image 
falling  In  several  pixels  would  be  less,  sometimes  twenty-five  percent  of  the 
signal  that  would  result  if  the  Image  had  fallen  in  one  pixel,  as  will  be 
shown  in  Section  II.  This  decreased  signal  is  the  reason  that  predicted 
satellite  detection  performance  and  false  alarm  rates  due  to  stars  would  be 
inaccurate  If  the  simplified  assumption  that  all  images  fall  in  one  pixel  Is 
made. 

Since  the  maximum  signal  Is  a  function  of  the  distribution  of  image 
photons  In  the  various  pixels;  the  number  of  pixels  that  an  image  occupies 
must  first  be  determined,  then  the  distribution  of  photons  In  those  pixels 


must  be  found.  The  maximum  signal  Is  due  to  the  pixel  with  the  most  photons, 
therefore  only  the  number  of  photons  In  this  pixel  have  to  be  known. 

The  energy  (number  of  Image  photons)  collected  by  a  pixel  can  be  deter¬ 
mined  geometrically  for  any  circular  Image  diameter  and  pixel  size  by  specify¬ 
ing  the  location  of  the  image  center.  If  the  Image  falls  randomly  on  the 
focal  plane  (no  tracking  of  Images),  the  probability  distribution  of  image 
centers  Is  uniform.  When  the  Image  is  confined  to  two  pixels,  the  area  that 
the  image  centers  can  occupy  is  limited  and  the  distribution  of  image  centers 
is  uniform.  The  expected  location  of  the  Image  center  is  therefore  the 
centroid  of  this  area  .  Once  the  expected  Image  center  location  Is  deter¬ 
mined  the  expected  energy  in  each  pixel  can  be  uniquely  determined.  For  the 
two  pixel  situation  the  expected  energy  in  each  pixel  is  50  percent  of  the 
total.  The  method  is  deceptively  simple  for  the  cases  of  the  Image  being  in 
one  and  two  pixels,  therefore  the  case  where  the  image  is  in  three  pixels  Is 
treated  in  detail  In  Section  II.  This  process  assumes  that  the  image  from  a 
point  source  Is  circular,  that  the  received  optical  signal  Is  uniform  over 
the  array  of  pixels,  and  the  pixels  are  of  equal  optical  responslvlty.  The 
method  described  In  this  paper  has  been  developed  for  the  case  of  square 
pixels,  no  dead  area  between  pixels,  and  for  uniform  photon  distribution  over 
the  circular  image  rather  than  Gaussian  or  Airy  functions.  The  presented 
method  can  be  extended  to  include  these  characteristics  at  the  cost  of  con¬ 
siderable  complexity  especially  for  the  photon  distribution  function.  Practi¬ 
cally  speaking  it  is  difficult  to  measure  the  photon  distribution  for  an 
optical  system  and  typical  measurements  are  of  the  nature  of  83  percent  of  the 
energy  in  a  seeing  disc  or  blur  circle  of  2  arc-seconds  diameter.  This  amount 


of  Information  describing  the  photon  distribution  function  is  appropriate 
for  the  presented  method. 


Since  these  Images  fall  at  random  locations  on  the  focal  plane,  the 
probability  that  each  Image  will  be  in  one,  two,  three,  etc.,  pixels  can  be 
calculated  based  on  the  geometrical  configurations  of  the  image  diameters, 
and  pixel  sizes.  Once  these  probabilities  are  known  they  can  be  multiplied 
by  the  expected  maximum  signal  for  one,  two,  etc.  pixels.  These  partial 
products  can  be  summed  to  give  the  expected  maximum  signal  due  to  any  image 
as  a  function  of  diameter  of  the  Image  and  pixel  size.  In  order  to  maximize 
the  SNR, the  design  of  most  electro-optical  systems  typically  results  in  the 
pixel  being  larger  than  the  expected  image  diameter  and  therefore  an  image 
rarely  falls  onto  more  than  four  pixels.  Due  to  poor  atmospheric  "seeing" 
conditions  the  image  diameter  will  grow  and  the  occurrence  of  the  image  in 
five  or  six  pixels  will  increase.  The  equations  required  to  determine  the 
expected  maximum  signal  and  the  probabilities  when  the  image  is  in  one  to  six 
pixels  are  given  in  the  next  section.  Plots  (Figs.  5-8)  are  given  that  allow 
the  expected  signal  to  be  determined  without  recalculating  the  equations 
presented. 

The  expected  maximum  signal  can  then  be  used  in  SNR  calculations  to 
predict  the  system  performance.  This  method  can  also  be  used  to  determine 
what  pixel  size  will  give  the  highest  SNR.  This  is  Illustrated  in  section  III. 
Although  the  electro-optical  detection  of  satellites  against  a  sky  background 
containing  stars  was  the  problem  that  motivated  this  work,  this  method  is 
applicable  to  other  detection  problems  where  circular  images  fall  on  square 
pixels. 


II.  EXPECTED  ENERGY  DERIVATION 

The  number  of  photo-electrons  due  to  a  point  source,  produced  by  an 
Imaging  device  at  the  focal  plane  of  a  telescope  is: 

I  =  p  a  q  t  (1) 

where:  p  =  photons  per  second,  per  unit  area,  from  a 
point  source, at  the  telescope  aperture, 
a  -  product  of  the  telescope  aperture  area  and 
the  telescope  optical  transmission, 
q  =  imaging  device  quantum  efficiency, 
t  =  integration  time. 

The  calculation  of  p  for  solar  spectrim  point  sources,  through  the  earth's 
atmosphere,  can  be  readily  accomplished  using  reference  (2)  by  R.  Weber. 

Equation  (1)  gives  the  total  number  of  photo-electrons  received,  but 
cannot  be  used  directly  to  calculate  the  SNR  because  the  calculations  depend 
upon  the  distribution  of  energy  in  x  number  of  pixels.  If  all  the  energy  is 
contained  In  one  pixel,  the  SNR  will  be  higher  than  if  the  energy  is  divided 
between  several  pixels. 

The  definition  of  SNR  is  the  peak  signal  divided  by  the  rms  noise.  The 
maximum  expected  energy,  as  described  in  Section  I,  is  the  peak  signal  due  to 
an  Image  falling  on  discrete  pixels.  The  threshold  value  used  for  detecting 
Images  (or  any  electrical  signal)  Is  determined  by  measuring  the  rms  noise  of 
the  Input  signal  without  an  Image  and  then  selecting  a  multiple  of  the  rms 
noise  as  the  threshold  value.  A  detection  system  can  be  characterized  by 
measuring  its  performance  vs.  the  SNR  of  the  input  signal.  The  SNR  of  the 


electro-optical  camera  is  therefore  an  Important  parameter.  This  section 
describes  the  maximum  energy  derivation  and  Section  III  demonstrates  how  this 
is  used  to  estimate  SNR  for  an  electro-optical  system. 

The  distribution  of  energy  from  a  point  source  among  the  Imager  pixels 
is  dependent  upon  the  size  of  the  circular  image,  the  size  of  the  pixels, 
and  the  centroid  of  the  image  on  the  focal  plane.  Assuming  no  optical  dif¬ 
fraction,  the  centroids  fall  randomly  on  the  focal  plane,  unless  the  images 
are  maintained  in  the  same  locations  by  tracking.  Since  it  is  equally  likely 
that  an  image  centroid  will  fall  on  any  location  within  the  array,  the  proba¬ 
bilistic  behavior  of  the  system  may  be  derived  via  geometric  analysis. 

The  well-known  expression  for  the  expected  value  of  a  discrete  random 
variable 

E(s)  =  E  S„  P„.  (2) 

may  be  interpreted  in  terms  of  maximum  values.  Then, 

Sn  =  the  maximum  signal  possible  from  a  single  pixel  within 
the  group  for  n  pixels  containing  all  of  the  signal. 

Pn  =  the  associated  probability  when  all  of  the  signal 
is  contained  in  those  n  pixels. 

E(s}  =  the  expected  maximum  signal. 

When  all  of  the  I  photo-electrons  are  contained  in  one  pixel,  Sj  s  I;  and  two 
pixels,  Sg  s  1/2. 

The  method  used  to  obtain  S3,  S^,  Sg,  and  Sg  is  identical,  therefore  only 
calculation  of  S3  wIV  be  ey  .ained  in  detail.  A  point  source  circular  image 
of  diameter  D  can  be  co. .ained  in  three  square  pixels  of  side  L  in  a  finite 
number  of  positions.  As  long  as  D  and  L  are  in  the  same  units  all  the 


following  equations  are  valid.  The  position  of  the  image  can  be  uniquely 


specified  by  the  location  of  the  image  center,  since  D  and  L  are  known. 

Figure  1  shows  the  possible  locations  of  the  image  center  when  an  image  is 
contained  in  any  three  pixel  combinations.  Since  the  image  center  can  only 
be  in  the  four  indicated  areas,  and  since  these  are  identical  and  symnetric, 
further  examinations  of  only  one  of  these  areas  is  sufficient.  Figure  2  shows 
one  of  these  areas  and  its  relationship  to  D  and  L.  Therefore  the  expected 
location  of  the  image  center  is  the  centroid  of  the  area  shown  in  Fig.  2 
(point  k).  Knowing  the  expected  location  of  the  image  center  and  the  image, 
the  ratios  of  area  in  the  three  pixels  to  the  total  image  area  can  be  calcu¬ 
lated.  These  ratios  multiplied  by  I  give  the  expected  signal  in  each  of  the 
three  pixels.  The  largest  of  these  ratios  is  the  only  one  of  interest  because 
if  it  crosses  a  detection  threshold,  the  point  source  is  detected.  If  the 
largest  signal  among  the  three  pixels  doesn't  cross  the  threshold,  the  other 
two  will  not,  and  there  will  be  no  detection. 

A.  Derivation  of  S3 

Due  to  geometric  contraints  there  are  only  three  ranges  of  D  and  L 
that  have  different  expressions  for  S3;  D  <  L,  L  <  D  _<  1.17L,  and  1.17L  <  D. 
Range  1.17L  <  D  results  in  S3  =  0  because  and  image  of  this  diameter  cannot 
fit  in  three  pixels.  Figure  3  shows  the  variables  used  to  calculate  S3,  x  is 

the  distance  from  the  pixel  center  to  the  expected  location  of  the  image 

center  (point  k).  At  is  the  area  of  the  image  and  A2  is  the  area  of  the  image 

not  in  the  pixel  with  the  maximum  energy.  For  the  range  0  <  D  <  L: 

A*  -  2A~ 

S3  =  - -  I  =  0.87  I  (3) 


At  i s  def i ned  the  same  for  all  Sn .  For  the  range  L  <  D  <  1 . 1 7L : 


where:  a  =  cos"^  (*~d/?~) 

x  -  (l/2  (L-D/2)2  -  USflfi  -  (f  -  f)  sincl)1/2 
z  =  L  -  D/2  -  x  cos  45° 


The  values  of  S3  for  this  range  vary  from  0.87  I  to  0.83  I. 

B.  Derivation  of  S4 

There  are  five  ranges  of  D  and  L  for  S4;  D  <  L,  L  <  D  <  1.17  L, 
1.17  L  <  D  <  1.25  L,  1.25  <  0  <  2L,  and  2L  <  D.  S4  =  0  for  2L  <  D.  For  the 
range  D  <  L: 

S4  =  .25  I  (5 

For  the  range  L  <  D  <  1.17L: 

S„  =  (.25  E,  +  E0)  I 


(6 


A2  '  (r-  ■  ?“)  cos  9  +  (6  -  s1n  w) 

.  l2  .  Di  .  (10^)  n  d2 

« ■  n--1  (H*) 

w  =  cos"1  (L/D) 
u  *  cos-1  ( . 7L/D) 

Ah2  =  At  -  (U^02W)  n  D2  +  (.35  sin  u  +  sin  w)  ^ 

\ 

The  values  of  for  this  range  vary  from  0.25  I  to  0.30  I.  For  the  range 
1.17  L  <  D  <  1.25  L: 


s4"  <-25Ei  E?>  1 

where:  A-j  =  (  L  -  ^)2 

A2  =  L2  -  (1/2  +  1/4  sin  w)  DL  -  (9Q--ygg  ~--a)  n  D2 
-  (L/2  -  D/4)  D  sin  a 
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cos 


.77L 

~r 


Ah?  ~  K  ~  n  D2  +  (.35  sin  v  +  sin  w)  2^=- 


The  values  of  S4  for  this  range  vary  from  0.25  I  to  0.28  I.  For  the  range 
1.25  L  <  D  <  2L: 


S4  =  .25  I  (8) 

C.  Derivation  of  Sg 

There  are  six  ranges  of  D  and  L  for  S^;  D  <  L,  L  <  D  _<  1 .17  L,  1.17  L 
<  D  <  1.25,  and  1.25  L  <  0  <  1.41  L,  1.41  L  <  D  <  2L  and  2L  <  D.  S5  =  0  for 
ranges  0  <  L  and  2L  <  D.  For  the  range  L  <  D  <  1.17  L: 

ss  ■  (if  Ei +  sf  £2) 1  <») 

where:  Ahl  =  -iny^o^--^  D2n  +  .75L  (D2  -  L2)1/2 

Ah2  ■  At  -  TO  "  °2  +  L  <d2  -  l2>'/2 

A,  -  Pmr*)  «  ♦  f- *'»  -  -  (f  -  £)  cos9 

ft2  *  (?  '  ?) 

The  values  of  Sg  for  this  range  vary  from  0.57  I  to  0.78  I.  For  the  range 
1.17  L  <  D  <  1.25  L: 


c  +  Ah2 

E'  *T 


(10) 
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where:  Aj  *  sin  w  +  (j  sin  a  -  L  +  £•]  /tan  a 

+  (gO-^zJwjn  D2  -  (L  -  |)2  tan  (90  -  a) 

A2  *  (^T")  '  (a720~)  n  D2  +  1.41  -  £  )  sin  (a  -  45) 

The  values  of  $g  for  this  range  vary  from  0.76  I  to  0.78  I.  For  the  range 
1.25  L  <  D  <  1.41  L 


-  /Ahl  c  +  Ah2  F  \ 

(*r  E>  *r  Ev 


where: 
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The  values  of  Sg  for  this  range  vary  from  0.73  I  to  0.80  I.  For  the  range 


1.41  L  <  D  <  21: 


■(N 


n2  n  ?  1/2 

where:  Aj  *  DL  -  \  -  (L  -  |)  (DL  -  L^) 


/ 180  -  2a  \  _  n2 
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The  values  of  Sg  for  this  range  vary  from  0.58  I  to  0.73  I. 

D.  Derivation  of  Sg 

There  are  five  ranges  of  D  and  L  for  Sg;  D  £  L,  L<  D  <  1.17  L,  1.17  L 
<  D  <  1.25  L,  1.25  L  <  D  <  2L,  2L  <  D.  Sg  *  0  for  ranges  D<  L  and  2L<  D. 
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For  the  range  L  <  D  <  1.17  L: 


M*;1 


where:  Ah  =  D2n  +  ^  (D2  -  L2)1/2 

The  values  of  Sg  for  this  range  vary  from  0.47  I  to  0.50  I. 

For  the  range  1.17  L  <  D  <  1.25  L: 

where:  Ah]  =  (^g*)D 2n  +  (D2  -  L2)V2 

Ah2 =  °2n  +  (L  - 1)2  +  (r  -  ?  )  s1n  a 

Ai  =  ih  °Zjl  -  k  (°2  ’  l^,Z 

a2  =  ((-w*)  d2ji  ■ 2-82  {LD "  sin  (45  ■  9))/  4L* 

j  -  tan-1  ((D2  -  L2)V2/l) 

The  values  of  Sg  for  this  range  vary  from  0.45  I  to  0.47  I. 

For  the  range  1.25  L  <  D  <  2L: 


where:  Ah  *  (L  -  |)  +(-4£-  -  £  )  sin  a  +  (^ygg-)  \ 

The  values  of  Sg  for  this  range  vary  from  0.28  I  to  0.45  I. 


The  method  used  to  derive  the  expressions  for  Pn  Is  similar  to  that  used 
to  obtain  expressions  for  Sn.  When  an  image  center  falls  at  random  on  the 
focal  plane  the  is  equal  to  the  area  that  the  center  can  occupy  while  the 

image  is  contained  in  exactly  n  pixels,  divided  by  the  area  of  one  pixel. 
Figure  1  shows  this  area  for  n  *  3.  The  expressions  for  Pn  are  a  function  of 
D  and  L,  and  change  for  various  ranges  of  D  and  L. 


P1  =  (L  -  D)2/L2  D  <  L 

Pj  =  0  L  <  D 

P2  =  2D  (L  -  D)/L2  D  <  L 

P2  =  0  L  <  D 

P3  =  D2  (1  -  H/4)/L2  D  <  L 

P3  =  (d2  (1  +  sin  a  -  JI  (^fgg))  -  DL  (4  +  2  sin  a)  +  4L2)/L2 

L  <  D  <  1.17  L 

P3  =  0  1.17  L  £  D 

P4  =  D2JI/4L2  D  <  L 

P4  =  ^DL  (6  +  2  cos  g  +  2  sin  a  -  sin  2  e)  -  4L2 

-  D2  (2  +  cos  g  +  sin  a  +  (45  -  2e  -  g  +  a))/L2 

L  <  D  <  1.17  L 

where:  g  =  sin  |(2L  -  D)/D^ 

P4  =  (d2  (1  +  sin  a  -  yjjg-  (90  -  2e  -  a)  -  DL 
(6  +  sin  2e  +  2  sin  a)  +  8L2  )  /L2 

1.17  L  <  D  <  1.25  L 


(16) 

(17) 

(18) 

(19) 

(20) 

(21) 
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(4L  -  4DL  +  D^/L1"  1.25  L  <  D  <  2L 

0  2L  <  D 

0  D  <  L 

(o2  (1  +  n/180  (90  -  a  -  4e)  +  cos  a)  + 

DL  (2  sin  2e  -  2  -  2  cos  a)  +  L2)/L2 

L  <  D  <  1.17  L 

(d2  (1  +  n/90  (67.5  -  a  -  2e)  -  1.41  sin  (a  -  45)) 
+  DL  (2  sin  (2e)  -  2  +  2.83  sin  (a  -  45))  +  L2 
-  (2L  -  D)2  tan  (90  -  a)  +  (D  sin  a  -  2L  +  D)2/ 


tan  a)/il 


1.17  L  <  D  <  1.25  L 


(DL  (4  -  sin  (2e))  +  D2  ( TK180  (2e  +  a  -  135)  -  2) 
+  L2  -  (4L  -  2D)  (DL  -  L2)1/2)  L2 


1.25  L  <  D  <  1.41  L 


=  (4DL  -  D2  (2  +  11/360  (180  -  2a))  -  (4L  -  2D) 

*  « 


(DL  -  L' 


1,1/2)/ 


1.41  L  <  D  <  2L 
2L  <  D 


D  <  L 


(tan'1  (2  (D2/4  -  L2/4)1/2/L) n  D2/180  -  2L 
(D2/4  -  L2/4)1/2)/L2 


L  <  D  <  1.17  L 


=  (d2  n  /360  (90  -  2g  +  2  tan  _1  (2  (D2/4  -  L2/4)1/2/L)) 
-  2L  (D2/4  -  L2/4)1/2  -  2.83  (DL  -  D2/2)  sin 
(45  -  g))/L2 

1.17  L  <  D  <  1.25  L 


p6  *  ({°  "  L/cos  9)2  sin  9  cos  g  -  L2  (4  +  tan  g) 

+  DL  (2  cos  g  +  2  sin  2  e)  +  D2  (1  +  n  /180) 
(45  -  2g  -  4e  -  cos  g)j/L2 


1.25  L  <  D  <  21 


Vo 


2L  <  D 


The  equations  for  Pn  where  n  =  1  to  6  can  be  used  to  calculate  the 


probability  of  an  Image  of  diameter  D  being  In  1  to  6  pixels  of  side  L  for 


any  combination  of  D  and  L.  Figure  4  Is  a  plot  of  Pn  where  L  =  4  and  D  varies 


from  0  to  8.  An  Image  D  *  3  falling  on  a  focal  plane  with  pixels  of  side 
L  *  4  has  the  probability  of  being  In  one  pixel  P]  3  0.06,  two  pixels  P?  *0.38, 
P3  *  0*12,  P^  =  0.44,  Pg  *  0,  and  Pg  *  0.  The  sum  of  the  probabilities  Is  also 
plotted  In  Fig.  4  for  0  *  0  to  1.25  L.  At  this  bound  the  probability  of  the 
image  being  In  seven  pixels  begins  to  be  greater  than  zero,  and  at  1.41  L  the 


probability  of  being  In  eight  pixels  begins  to  be  greater  than  zero. 


Review  of  equations  (3)  through  (15)  shows  that  Sn  can  be  rewritten  as: 


S  *  k  I 
n  n 


This  allows  equation  (2)  to  be  rewritten: 


E(s}  »  I  k  I  P  *  I  E  k  P 
n  n  n  n  n  n 


From  equation  (31)  a  normalized  expected  maximum  signal  can  be  expressed  which 


shows  the  dependence  on  D  and  L,  and  Independence  on  I. 


NOR  E(s>  =  E{s  VI  *  I  tc  P„ 

_  n  n 
n 


ru 


Fig.  4.  Probability  of  an  Image  falling  into  n  pixels  for  image 
diameters  (D)  0  to  8  and  for  pixel  size  L=4. 


The  values  of  equation  (32_)  can  be  viewed  as  attenuation  factors  of  I  that 
result  In  the  expected  signal  E(s>.  Equation  (32)  Is  plotted  in  Fig.  5  to  8 
for  many  values  of  D  while  L  varies  from  D/1.25  to  15.  The  lower  limit  of  L 
is  again  due  to  the  image  not  being  considered  in  more  than  six  pixels.  The 
values  of  NOR  E{s)  increase  as  L  increases  which  by  itself  suggests  that  big 
pixels  give  more  signal.  However  while  E{s>  is  the  numerator  of  the  SNR,  the 
denominator  Increases  at  a  much  faster  rate  and  results  in  the  important  SNR 
decreasing  with  large  pixel  size.  This  is  discussed  in  the  next  section. 

The  relative  Increase  in  NOR  E{s)  when  L  <  D  is  due  to  Pg,  Pg,  Sg  and  Sg 
increasing.  Although  Pg  and  Pg  would  be  expected  to  Increase  as  D  becomes 
larger  than  L  It  Is  not  intuitive  that  Sg  and  Sg  should  Increase.  They 
increase  due  to  the  geometric  contraints  Imposed  when  a  circular  energy  dis¬ 
tribution  is  in  5  or  6  square  pixels.  This  causes  the  "tails"  on  the  curves 
in  Figs.  5  to  8.  Values  of  Sg  range  from  0.73  I  to  0.80  I  and  Sg  from  0.45  I 
to  0.50  I.  Intuitively  several  people  have  suggested  that  Sg  should  be  1/5  I 
and  Sg  should  be  1/6  I.  A  simple  graphical  check  will  show  that  this  is  not 
possible. 

Using  the  values  of  NOR  E{s}  plotted  In  Figs.  5  to  8,  the  E(s>  may  be 
calculated  for  any  I  without  having  to  calculate  equations  3  through  29  for  a 
particular  combination  of  D  and  L.  D  and  L  must  be  in  the  same  units. 


III.  USING  THE  EXPECTED  ENERGY 


The  maximum  expected  energy  from  a  point  source  can  now  be  determined,  and 
used  to  calculate  SNR. 

E{  s } 

SNR  =  TT7I  (33) 

(B  +  N2) 

where:  B  =  mean  photo-electrons  due  to  background 

N  =  RMS  noise  electrons  of  imaging  device 
The  calculation  of  B  is  similar  to  that  of  I. 

B  =  b  a  q  t  L2  (34) 

where:  b  =  photons  per  second,  per  unit  area,  per  angular 

area  of  sky  at  the  telescope  aperture. 

The  calculation  of  b  can  be  accomplished  using  reference  (2).  The  numerator  of 
equation  33  increases  as  L  increases,  as  shown  in  Section  II  using  the  normal¬ 
ized  expected  signal.  The  denominator  of  equation  33  also  increases  with  L 
and  to  demonstrate  this  increase  and  dependence  of  L  the  denominator  divided 
by  (b  a  q  t  +  N2)^2  is  plotted  in  Fig.  9  for  b  =  792,  a  =  0.46,  q  =  0.4, 
t  =  .6,  and  N  =  20.  This  results  in  a  normalized  plot  of  the  SNR  denominator, 
although  not  normalized  to  one. 

The  input  SNR  required  to  produce  a  specified  performance  of  detection 
(probability  of  detection  and  false  alarm  rate)  is  a  useful  way  in  which  to 
evaluate  the  performance  of  an  electro-optical  system.  In  essence,  the  system 
is  viewed  as  a  "black  box".  The  performance  of  a  detection  system  as  a  "black 

box"  can  be  measured  as  a  function  of  the  input  SNR.  In  many  situations, 

4  5 

a  minimum  acceptable  SNR  is  from  6  to  8.  An  automatic  detection  system 


1 

* 

! 

f 

> 

i 

can  detect  targets  with  a  probability  of  detection  of  .95,  a  false  alarm  rate 
I  -3 

1  of  4.8  *  10  ,  when  the  SNR  is  as  low  as  3. 

Plotting  SNR  using  equations  (1)  through  (34)  shows  the  dependence  of 
SNR  on  L  for  several  values  of  D.  Figures  10  through  18  are  plotted  with 
j  typical  parameter  values  which  are  listed  on  each  plot.  Area  is  in  meters 

;  squared,  time  is  in  seconds,  0  and  L  are  in  arc-seconds  and  all  other  units 

are  specified  for  equations  (1)  and  (34).  All  plots  start  at  L  *  D/1.25 
where  the  image  is  in  six  pixels  or  less.  The  tails  on  the  plots  when  L  <  D 

I 

;  is  again  due  to  the  geometric  contraints  imposed  when  a  circular  image  is  in 

i 

;  5  or  6  pixels  as  discussed  in  Section  II. 

Figures  10  through  14  are  plotted  with  p  (photons  from  a  point  source) 
and  b  (photons  from  the  background)  changing,  while  the  parameters  a,  q,  t, 
and  N  are  fixed.  This  is  the  case  in  an  operating  electro-optical  system 
where  the  satellite  and  sky  are  changing  in  intensity  over  time.  Figures  15 
through  18  have  all  parameters  fixed  to  match  Fig.  10  except  for  one  parameter 
in  each  figure;  i.e.,  a  in  Fig.  15,  q  in  Fig.  16,  t  in  Fig.  17,  and  N  in  Fig. 
18.  This  demonstrates  the  degree  that  the  electro-optical  system  designer 
can  affect  SNR. 

'ii 

Although  the  shape  of  the  SNR  plots  in  Figs.  10  through  18  change,  depen¬ 
dent  on  the  parameter  values  used,  it  is  clear  the  there  Is  an  optimal  pixel 
size  to  maximize  SNR.  The  importance  of  pixel  size  Is  frequently  overlooked 
while  the  other  controlable  parameters  a,  q,  and  t,  are  maximized  and  N  mini¬ 
mized.  There  is  an  intuitive  feeling  that  the  pixel  size  should  be  large  In 
the  hope  that  a  point  source  Image  will  fall  in  one  pixel  and  therefore  give 
a  large  signal.  While  at  the  same  time  it  is  clear  that  a  large  pixel  will 
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Fig.  12.  SNR  vs.  L,  b  =  4159. 
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Fig.  17.  SNR  vs.  L,  t  =  .50. 
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collect  more  background  photons  which  will  increase  the  noise.  The  method  of 
estimating  the  expected  maximum  signal  derived  in  this  paper  and  the  demon¬ 
stration  of  incorporating  this  in  SNR  calculations  should  enable  the  selection 
of  pixel  size  to  be  optimal. 

The  system  designer  may  have  some  control  over  D  (image  diameter)  by 
telescope  specification  and  by  site  selection,  however,  usually  the  atmosphere 
results  in  a  range  of  D.  Figures  10  through  18  demonstrate  that  it  may  be 
better  to  maximize  the  SNR  for  small  D  by  selecting  a  small  pixel  size,  and 
then  be  penalized  when  the  "seeing"  is  bad  (D  increases),  rather  than  select¬ 
ing  a  large  pixel  size  and  always  have  low  SNR. 

Figures  10  through  18  give  a  concrete  illustration  of  how  the  correct 
pixel  size  can  be  selected,  what  SNR  will  result,  and  if  the  optimal  pixel 
size  cannot  be  selected  (due  to  practical  constraints)  what  penalty  will  be 
paid  in  SNR.  In  addition,  it  will  be  evident  that  the  electro-optical  system 
will  provide  sufficient  or  insufficient  SNR  for  the  detection  system,  human 
or  electronic. 
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